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I. INTRODUCTION 

^ ' Description of weak interaction processes in nuclei is often a challenge for nuclear structure models. Numerous 

Ph '• calculations of the nuclear double beta {/3(3) decay amplitudes well illustrate this statement (see, e.g., Refs. [1]). 

Uncertainties in theoretical calculations of the Gamow- Teller (GT) two-neutrino /3/3 (2z^/3/3) decay amplitude Af^^ 

fvq , have stimulated experimental studies of the GT^^^ strengths for the transitions to the intermediate 1+ states virtually 

^ ' excited in the decay process (see, e.g. Refs. 0-0 )• 

■TJ" I As a double charge-exchange process, 2i/(Oi^)/3/3-decay is enhanced by nucleon pairing which is due to the singlet 

r ■ ■ part of the particle-particle (p-p) interaction. The discrete quasiboson version of the quasiparticle RPA (pn-dQRPA) 
^^ I which accounts for the nucleon pairing is usually applied to calculate the /3/3-decay amplitudes in open-shell nuclei 
. . [l|. In spite of differences in the model parametrization of the nuclear mean field and the residual interaction in 
l/~) ' the particle-hole (p-h) and p-p channels, all pn-dQRPA calculations reveal marked sensitivity of the amplitude Mq^j, 
^^ to the ratio gpp of the triplet to siriglet p-p interaction strengths. Physical reasons for such a general feature of all 
calculations were analyzed in Ref. [5|, and the sensitivity was shown to reflect violation of the Wigner spin-isospin 
SU(4) symmetry in nuclei. An identity transformation of the amplitude into the sum of two terms was used in 
Ref. J5|. One term, which is due to p-p interaction only, depends linearly on gpp and vanishes at gpp = 1 where the 
SU(4)-symmetry is restored in the p-p sector of a model Hamiltonian. The second term is a smoother function of gpp 
at gpp ~ 1, but exhibits a quadratic dependence on the strength of the mean- field spin-orbit term, which is the main 
d ' source of breaking of the spin-isospin SU(4)-symmetry in nuclei. 

Understanding of general properties of the amplitude Mq^j, helps to improve the reliability of evaluation of the 
/3/3-decay amplitudes. For a quantitative analysis, we use here an isospin-selfconsistent pn-continuum-QRPA (pn- 
cQRPA) approach, which was initially proposed in Ref. [yj to describe different strength distributions in single-open- 
shell nuclei. In the reference the full basis of the single-particle (s-p) states was used in the p-h channel along with 
the Landau-Migdal forces, while the nucleon pairing was described within the simplest version of the BCS-model 
based on a discrete s-p basis. A rather old version of the phenomenological isoscalar nuclear mean field (including 
the spin-orbit term) was used in Refs. [a, |y], as well. A further development of the pn-cQRPA approach applied to 
description of the /3/3-decay matrix elements has been given in Ref. Q . In this reference realistic (zero-range) forces 
have also been used in the p-p channel to describe the nucleon pairing within the BCS model realized on a rather 
large discrete-|-quasidiscrete s-p basis. 

In the present work along with the detailed formulation of the pn-cQRPA approach we use a modern version of the 
phenomenological isoscalar mean field (including the spin-orbit term) deduced in Ref. [8| from the isospin-selfconsistent 
analysis of experimental single-particle spectra in double-closed-shell nuclei. The last version of the approach and 
some applications to the GT observables are briefly described in Ref. [^. 

The paper organized as follows. In Sect. 2 starting from the coordinate representation of the pn-dQRPA equations, 
we formulate a pn-cQRPA approach. The expressions for different type strength functions are also given. A realistic 
model Hamiltonian and the calculation scheme are described in Sect. 3. In Sect. 4 the results of pn-cQRPA 
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calculations of the GT^ strength functions and Af^^ amplitudes for the number of nuclei are presented and compared 
with the available experimental data. Conclusions are drawn in Sect. 5. 

II. VERSIONS OF THE PN-QRPA: DISCRETE AND CONTINUUM 

A. pn-dQRPA equations in the coordinate representation 

The first step towards formulation of a pn-cQRPA approach is transformation of the pn-dQRPA equations to the 
coordinate representation (see also Refs. [y, 0, Q)- The energies w^ and wave functions \s^J^M) of the states in 
the isobaric (odd-odd) nucleus are usually obtained within the quasiboson version of the pn-dQRPA as a solution 
of a system of homogeneous equations for the forward and backward amplitudes Xf^^ and Y/^^ related to /3^ 
charge- exchange excitations of an even-even parent nucleus. Here, J and M are the total angular momentum and its 
projection and tt = (—1)'^ is the parity of the intermediate states (L = 0, 1, 2 . . . and S" = 0, 1 are the transferred 
orbital momentum and spin, respectively). 

The pn-dQRPA eigenenergy solutions w^ are related to the excitation energies EifJ measured from the ground-state 
energy Eo{Z ± 1, TV =F 1) of the corresponding daughter nuclei as: 

C.W = u;, ± (mp - /in) = i^i^s^ + QI^^ (1) 

Here, /ip(„) is the chemical potential for the proton (neutron) subsystem found from the known BCS equations, 

Qj^ = £b{Z,N) — £b{Z±l,N^l) are the total binding-energy differences, cus are the excitation energies measured 

from Eo{Z, N) — J2a "^a ~ -~£b{Z, N) {rria is the nucleon mass). The energies tOs are usually described by a model 
Hamiltonian. 

For definiteness, we consider further excitations in the /3~ channel. The amplitudes X^^^ and Yf^^ are based on 
the discrete s-p levels n and z/ for protons and neutrons, respectively (each set tt and u contains the s.p. quantum 
numbers rt,., I, j). The excitations in the /?+-channel are described within the pn-dQRPA by the same equations with 
substitution p ^ n {tt -n- u). 

Following Ref. [y] , we transform the system of equations for the amplitudes X and Y into the coordinate represen- 
tation. Transformation is performed in terms of the four-component radial transition density r''^ gj^^ (r) defined as 
follows: 

P^f '(r) - \gi':^{r)TjLSM{n), gii'{r) ^ ^ g^^^^^ xit'{r), (2) 



J,. 







^TT^U^TTU l~ ^7T ^V^^TTV 



where u, v are the coefficients of Bogolyubov transformation; xH^i^) = ^mIu) X7r(^)Xi'(^) '^ith \/2 J + 1 t-Jh^^. — 
{{tt)\\Tjls\\{'^)) being the reduced matrix element of the spin-angular tensor TjLSM{n = r/r); (tt) — Jt^I-k 
ii.^) = JutIu) and r~^XTr{T') {f~^Xv{i'))i being the s-p proton (neutron) quantum numbers and radial wave func- 
tions, respectively. Hereafter we shall sometimes omit the superscript "JLS^^ and/or the subscript "s" when it does 
not lead to a confusion. 

According to the definition ([2]), the elements gi, g2, gs, gi can be called the particle-hole, hole-particle, particle- 
particle and hole-hole components of the transition density, respectively. In particular, the transition matrix element 
to a state \s, J^M) corresponding to a probing particle- hole operator 

VJlIm =yyL {ra)TjLSM (Ha) T^") (3) 



is determined by the element gi as: 



Assuming the residual interaction F^ in the p-h (K = 1,2) and p-p {K — 3, 4) channels to be momentum-independent, 
one can represent the system of homogeneous pn-dQRPA equations (primarily written for the amplitudes X and Y) 



{s, rM\VJll,,\0) = I e(^^(r)yL (r) dr. (4) 



in the equivalent form: 

K' •' 

Here, {rr')~^F^{r,r') is the radial part of the residual interaction, the 4x4 matrix (rr')~^j4/f/f/(r, r', w) is the 
radial part of the free two-quasiparticle propagator hereafter taken diagonal on L, S (the so-called "symmetric" 
approximation @): 

AKK'{ri,r2,Uj) =J2 XTri'{ri)X7n^{r2) AKK',7r,^{(^), AKK',Tru ^ Ak'K,ttv, (6) 

ulvl ulvl ulul vlvl 

All,TTU — -^ —^ , ^33, TT 
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A\2,-nv = -^34,7ri/ = U^V.„Vt,U,, 



Ld — Et^i, LO + E-,^1, J 



^13,7r,y = UuV^{ 2_ 1 "" ), Ai4_^^ = -U^U^(- 



A22,wiyiuj) = ylii^Triy(-w), Aaa,t,u{i^) = ^33,7ri/("W) , A237riy(w) = Ai4,7riy (-1^), A24,7riy (i:^) = Ais^jri/ (-w) , 

with E'ttz^ = ^'tt + E^, where Et^{E,^) is the proton (neutron) quasiparticle energy. Excitations in the /3+-channel 
are described within the pn-dQRPA by Eqs. (O,® with the substitution p <->■ n (tt o z/). The relationship 
Akk' .utt{'^) — Akk' .i:u{~^)i which follows from Eq. (jG)), allows one to reduce pn-QRPA description of excitations 
in the /3+-channel to pn-QRPA description of excitations in the /^"-channel. In particular, because of equalities 
^ii,i/7r(w) = Aii.Tr,y(— w) = A22,7ri/(w), the transition matrix element to the state |s, J^ M) corresponding to a probing 
operator Vj^^^j (given by Eq. ([3]) with substitution t^ ' -> t^^') is determined by the transition-density element 

Q2'- {^■•J^^Wjlsm\^) ~ S (^2^s^{''')^L{f')dr. In such a description the operator Vj^^ji/ can be considered as the 
hole-particle one. 

B. QRPA effective operators and strength functions 

The next step towards taking the s-p continuum into account is consideration of the effective two-quasiparticle 
propagator (two-quasiparticle Green function) AKK'ir^r' ,u}), which satisfies a Bethe-Salpeter-type equation |7j|: 

Ai^i,{ry,.j) = Ai^i,{ry,uj)+Y.hKK''yr,,^)FUrur2^ (7) 

K" •' 

Here, the current energy uj is related to the excitation energy of daughter nuclei: ui^^' = w ± (/Xp — /i„) (compare 
with Eq.(l)). The spectral decomposition of the effective propagator contains all the necessary information about the 
pn-QRPA solutions: 

A,, in,r2,.)^^ '^-^,+,0 -L y + ^^^,0 ' ^^) 

s s 

Aii'^in,r2,Lj) = iff^(ri,r2,-c.), (9) 

A,, in,r2,.)^^ '^-^,+,0 -L ^ + ^,_,o - (10) 

s s 

In practical applications of the pn-QRPA method it is more convenient to use four-component effective fields 
^x[/](''>w) corresponding to the probing operator V{r) acting in the particle-hole (/ = 1) or hole-particle (/ = 2) 
channels. The effective fields are defined as follows: 



/ AKiir, r', w) V{r') dr'^Y. [ ^^^' (''' ^'^^Wk'[i] i^' , ^) dr' . 



(11) 



The system of equations for the effective operators follows from Eqs. (I7)). (|lip : 

yif,f (r, c) = Vdr)5Ki + J2 [ ^Kir, r') A^^^ (/, r", c.)^^,^!^ (r", a;) dr'dr". (12) 

K' •' 

Similar equations are given in Ref. [101: where they are obtained by the methods of the Migdal's finite Fermi-system 
theory. 

By definition, the strength functions corresponding to the probing operators of Eq. ^ arc: 

4LsH-E|(^-^"ll^iLsl|0)f'5(..-c.,). (13) 

s 

Making use of the spectral expansion Eq. Q, one can express the strength functions in terms of the two-quasiparticle 
effective propagator, or in equivalent terms of the corresponding effective fields: 

4lsH = -i(2J+ l)Im/ Vi(ri)iff^(ri,r2,L^)F(r2)dridr2 = 
= -i(2J+l)ImE/^L(r-i)A5'F(^i,'^2,^)l>ifif(r2,w)dridr2. (14) 

K 

= -ilm(2J+l)E/^L(^i)^2^K^(^i,^2,^)V'if2f(r2,w)dridr2. (15) 

K 

Within the pn-QRPA description of the /3/3-decay amplitudes it is useful to consider a "non-diagonal" strength 
function Q: 

s^iJ{u)^Y.^Q^\vy^l\\s,r){s,r\\v!,i>sm6{^-^s), (16) 

s 

where |0') is the ground state wave function of the final nucleus. Identifying the pn-QRPA vacuum |0') with that for 
|0), one gets the following expressions: 

4ls^M = -i(2J+l)Im/V^L(ri)4f ^)(ri,r2,L^)FL(r2)dridr2 = 

= ~^{2J+l)lmY,jVL{ri)A^^j^'^\r^,r2,uj)V^f,nr2,u:)dr,dr2. (17) 

K 

Bearing in mind applications of the method to describing the GT strength functions and 2j//3/3-decay amplitude 
Mgy, in this paper we limit concrete consideration by the case (JLS) = (101) and use for this case the probing 

GT operators Gj,j- = ^^ aMia)Ta ■ Here aM are the Pauli spherical matrices, the radial part of the corresponding 

external field of Eq.Q is VL=o{r) = s/A'k. The integrated strength functions Sqj, satisfy the Ikeda sum rule equal 

3(iV-Z). 

C. 2!//3/3-decay amplitude 

The nuclear GT~ amplitude for the 2i//3/3-decay into the ground state |0') of the final nucleus (Z-|-2,N-2) is given 
by the expression 

where oj^ = E'^ - \{Eq + Eq,) = E^^s + \{Q'i^ + Q'-^-'b)- To calculate M'^rj. within the pn-QRPA, the vacua |0) 
and |0') should be identified. As a result of such identification, one has tD^ — \{i^s +'^'^''i) ~ ^s in accordance 



with Eq.(l). Therefore, the amphtude ^TE\\ can be directly expressed in terms of the "non-diagonal" GT strength 
function of Eqs. (ITBl). ([T7|: 

M^^j,^ f uj-^sl7jr\uj)duj. (19) 

An alternative expression for Mq^j. is obtained in terms of the "non-diagonal" static polarizibility [3] : 

A^-y = -SV^^I^rC^^',^ = 0)V;[,-^(/,L. = 0)drdr'. (20) 



As mentioned in Introduction , the amplitude of Eq. ()18p can be decomposed in a model-independent way into two 
terms [5i]: 

M'aT = M^^ + M^^, M-^^ ^^i^^^ , (21) 

'^GTR 

where Ojgtr is the energy of the GT^ giant resonance (GTR), and EWSRqj, — ^ loSqj, {uj)duj is the "non- 
diagonal" energy-weighted sum rule. This sum rule, which is determined by the p-p interaction only, is mainly 
responsible for the sensitivity of the amplitude M^.^ to the ratio gpp of the triplet to singlet strength of the p-p 
interaction. The sum rule EWSRQrp depends linearly on gpp and vanishes at gpp = 1 where the SU(4) symmetry 
is restored in the p-p sector of a model Hamiltonian. The term Mqj, in Eq. (|21l) is a smoother function of gpp at 
gpp ^ 1, but exhibits a quadratic dependence of the strength of the mean- field spin-orbit term, which is the main 
source of violation of the spin-isospin SU(4)-symmetry in nuclei. 

D. Exact treatment of the single-particle continuum in the particle-hole channel 

The above given coordinate representation of the pn-dQRPA equations is fully equivalent to the version usually 
formulated in terms of the X and Y amplitudes, provided a discrete basis of s-p states is used. However, the use 
of the coordinate representation allows us to exactly take into account the s-p continuum in the p-h channel and, 
therefore, to formulate a version of the continuum-pn-QRPA (pn-cQRPA) (see also Refs. [6|,0,Im|)- Within this version 
the pairing problem is solved on a basis of bound and quasibound proton and neutron s-p states. To take the s-p 
continuum into account , the following transformations in Eq. ^ is performed: (i) the Bogolyubov coefficients v\, 
u\ and the quasiparticle energies E\ are approximated by their non-pairing values v\=0, u\ = 1 and E\ = e\ — fj, 
for those s-p states, which lie far above the chemical potential (i.e., {e\ — n) ^ A); (ii) the radial s-p Green function 
g{\)ir,r' ,e) — J2e (^ ~ ^ a + *0)^^xa('')xa('"') is used to explicitly perform the sum over s-p states in the continuum 
(A = Ea, (A); (A) = jA, l\)- As a result, we get from Eq. (|6]) the following representation of the properly transformed 
free two-quasiparticle propagator: 



1 



A I \ Y^ ^'^yW^X7r^/(^l)X7^Iy(7'2) , v;^ ,2 i \ i \ 2 I { , t^i \ 

Aii{ri,r2,uj) = 2^ + 2^ i(^)(^)Xi'('-i)x>.('-2) w^5(^)(r-i,r2,Ai^+w-£;„) 

-|- {p o n, TT <->■ (y, a; O — w} , 

Axi{r-Y,r2.,<^) = V" u^Vt,UyrVyrX-^i^{ri)x-^u{r2){ ; 

? ? 

Ai3{ri,r2,uj) = V X,n.(r-i)x7r,.(?'2)wi.wi.( ^ h — r^j 

+ X! ^l^)(u)Xu{ri)Xv{r2)uuV„g[^){ri,r2,uj + ^.p - E,,), 



2 2 

Ai4,{ri,r2,Uj) = - V X7ri/(?-l)X7n.(?'2)w^W7r( ^ 1 rTT-) 

7r,z-'< 

7r<,(i/) 

A33[ri,r2,UJ) = > X^>^(^i)X7ri'(r2)( — ^^ 



+ X] *W(>^)X7r(ri)X7r(r2)u^5(,.)(^l,^2,W--B^+/^n) 
7r<,(i') 

+ Yl i'{7r)iu)XAri)x.^ir2)ulg[^^iri,r2,uj-E^ + fip). (22) 

i^<,(7r) 

Here, tt < (i^ <) means tt < TTmax [v < Vmax)-, where TT^ax {vraax) is taken so that i?7r,„„^ = £7r,„„, - Mp (-£^i^™„x = 
Si^mox ~ Mn) with an acceptable accuracy, and the projected radial s-p Green function is: 

/ / N / N Y^ Xa(?'i)Xa(»-2) ,„„s 

, ^ ^ A 

The first sum in this expression is running over all the s-p states and can be calculated via a product of the regular 
and irregular solutions of the corresponding s-p Schrodinger equation. Such a representation of 5(a) (^i, ^2, e) has 
been used in Ref. [11] to formulate the very first version of the continuum-RPA. Other AKK'iTTf' ^ijj) components in 
Eq. ([22]) are obtained in the same way, as it is done in Eq. ([6]). Thus, Eq. (|22|) together with Eqs. (f23| and (fT2|) 
realize a version of the pn-cQRPA method. As applied to description of the Fermi and GT strength distributions in 
single-open-shell nuclei, the method was used in Ref. \g\. This pn-cQRPA approach has been applied to description 
of the /3/3-decay matrix elements in Ref. [Tj] . 

Within the pn-cQRPA the spectral decomposition of Eqs. ([5|)- (ITU|) remains valid, provided the continuum-states \s) 
(normalized to the (5-function of the energy) are taken into consideration. These states are characterized, in particular, 
by the particle-hole (hole-particle) energy-dependent transition densities which determine the strength functions as 
follows: 



S^j%{lo)^{2J+1) 



2 



VL{r)Qyi!sirMdr 



(24) 



These transition densities can be expressed in terms of the effective fields corresponding to the probing operators of 
Eq. (|3]) and their partners. Supposing the p-h interaction is of zero range (with the intensities F^ = F^ = 2F^'), one 
gets from Eqs. (fT2|) - ([T5|) the expressions: 

(_ 1 i2J+l)^^'vii[^'\r,.) (,, 1 (2J+l)V^y<^f^)(r,.) 

The intensities F may be radial-dependent. 

III. CALCULATION OF THE GT STRENGTH FUNCTIONS AND THE 2iy 13 fi-HEC AY AMPLITUDES 

WITHIN THE PN-CQRPA 

A. Model Hamiltonian: interactions in the particle-hole and particle-particle channels 

For quantitative realization of the relationships given in the preceding Section, we specify the nuclear Hamiltonian. 
We adopt here the realistic model Hamiltonian, which is an extended version of that used in Refs. [0, 0|- The 
Hamiltonian consists of the mean field and residual interactions. The partially selfconsistent mean field is described 



in the Subsect. B. Here, we specify interactions in the p-h and p-p channels, Hp-h and Hp^p, respectively. As Hp^h, 
we use the Landau-Migdal forces in the the charge-exchange channels: 

Hp^h = 2^(F° + F'aa(7i,)Ti-hl+^S{ra - n) + h.c, (26) 

a>b 

where the intensities F of the non-spin-flip (5* = 0) and spin-flip {S = 1) parts of this interaction are the phenomeno- 
logical parameters. We choose the zero-range p-p interaction in both the neutral (pairing) and charge- exchange 
channels: 



p—p 



1 Y: GHXX'X,X[Us{pL'f)'pi%, (27) 



JLSM 



where G^ are the intensities of the p-p interaction in the S* = 0, 1 channels, and (differently from the below-given 
equations) the sums over s-p states include the neutron and proton states, {XX'XiX'i)jls — /Xa a'^('')Xa^a' ('')75- 
The annihilation operator P'"^^ for a nucleon pair is given, as follows: 

P'lx' = Y. {JM\jxm,jx'm')ay^'ax^, (28) 

m,m' 

where a\„i {a\^) is the annihilation (creation) operator of a nucleon in the state with the quantum numbers Xm 
(rn is the projection of the particle angular momentum), {JM\j\m,j\im') is the Clebsh-Gordan coefficient. The p-h 
interaction of Eq. (I26p can also be expressed in terms of the operators a^ and a: 

Hp^H^2 J2 F''in,.n',.')jLs{Ait'yAi^^,+ (n^X (29) 



JLSM 
■ku-k' y' 



where the creation p-h operator A^ is given as follows: 

(Ai^)" ^ Y, {JM\j^m^,j,~m,}{~iy"-^"ai^^a,^^. (30) 

In solution of the pairing problem we simplify the p-p interaction (j27p for the 0"*" neutral channel, using the "diagonal" 
approximation: A = A', Ai = X[. The nucleon pairing is described with the use of the Bogolyubov transformation in 
terms of the quasiparticle creation (annihilation) operators a\„^ , aA777 (see, e.g., Ref. (iJl)- As a result, the following 
model Hamiltonian is obtained to describe charge-exchange excitations within the pn-QRPA: 

H ^ Ha^n + Ho,p + fJ-nN + fipZ + Hp-h + Hp-p, Hq = Y^Exa\,^axm- (31) 

Ar77 

Here, Ex — \/{ex — /i)^ -I- A^ is the quasiparticle energy. The chemical potentials /i and the energy-gap parameters 
Ax are determined by solving the BCS-model equations for even and odd subsystems (with the number of nucleons 
N{even) and N{odd) nucleons, respectively): 

N{even)^Y.^23x + l)vl, AA(e«en) = G" ^(2jv + l)(AAA'A')ooo^, (32) 



N{odd) = Y, {23X + 1)1^1 -t- {23Xo - IX + 1, (33) 

X^Xo 

Axiodd) = G°\ V (2jv + l)(AAA'A')ooo:^ + (2jao - l)(AAAoAo)ooo:^ j- (34) 



Here, Aq labels the level occupied by the odd qausiparticle. Following the standard procedure, we introduce different 
strengths G^ and Gp of the p-p interaction in the 0"'' pairing channels to reproduce the experimental pairing energies. 
The latter are expressed via the ground-state binding energies calculated within the BCS model: 

- Eb{even) = ^(2ja + l)vl{^i - Ex), (35) 

A 

- £i,{odd) = J2 (2JA + 1)«^(m - Ex) + (2jAo - 1)«^(a* - Ex,) + sxo- (36) 

The level Aq is chosen from the minimum of the calculated — f b, and/or accordingly to the quantum numbers of the 
odd-nucleus ground state. To roughly correct for double counting in evaluation of the total energy, we properly modify 
Eqs. ([55]) . (1551) : (^ — Ex) — J> {fi — Ex — ^Uxx)', £Xo ~^ ^Xo ~ ^UxgXo- Here Uxx is the diagonal matrix element of the 
nuclear mean field (Subsect. B.). 

The explicit expressions for the p-h and p-p interactions in terms of the quasiparticle (pn)-pair creation and 
annihilation operators are: 

iJp_, = 2 Y. F'i^'^^''^')jLs(Qi'jyQi^S + (p^n,n^A (37) 

JLSM ^ '^ 

4-P = - E G^(^^^V').Ls(^.r)V;^,t^, (38) 

JLSM 

7TU7T v' 

where 

A'J' - E {JM\j.m,j,m')a,„,,a.,n, Kf = {-!)'-'' Ai-'' . (40) 

mrn' 

The forward and backward amplitudes X and Y are defined by the corresponding matrix elements X^^^ = 

{s,J^\{A:l^^)^\0) and Y/J;'^ ~ {s,J'"\Ai^\0). The system of equations for these amplitudes can be explicitly de- 
rived using the model Hamiltonian of Eqs. ([3T|). ([37)) - (l40l) . Having transformed this system to the form given by 
Eqs. ©J®,®, one can relate the interaction F^{r,r'), entering the basic pn-QRPA Eqs. ®, p^ . to the model- 
Hamiltonian zero-range interactions of Eqs. (PTjl - lHUj) : 

F^{ry)=F^^±^, F|=i,2=2f^ fL^a = -G'. (41) 

rr 

Thus, within the model under consideration there are four phenomcnological interaction strength parameters F , G^ 
{S — 0,1) which will be specified below. 

B. Model Hamiltonian: the mean field and selconsistency conditions 

The nuclear mean field consists of the isoscalar (central and spin-orbit), isovector and Coulomb parts: 

U{x)^U„ix) + U,{x) + Uc{x), Uo{x) = U„ir) + Uso{r)lrT, (42) 

Ui{x) = ^v{r)T^'\ Ucix) = ^il-T<^'^)Uc{r), (43) 

where v{r) is the symmetry potential. Within the model the mean-field isoscalar part is treated phenomenologically 
(that seems well justified for description of isovector excitations), while the isovector and Coulomb parts are calculated 
selfconsistently. 



The realistic phenomenological isoscalar mean field is parametrized as follows: 

Uoir)=-Uofws{r,R,a), Uso{r) ^ Uso-^^, (44) 

r dr 

where fiys — {l+exp{^^^))~^ is the Woods-Saxon function with the radius R — tqA^^^ and the diffuseness parameter 
a. To avoid spurious violation of the isospin symmetry caused by the mean- field isovector part, the latter is calculated 
using the isospin-selfconsistency condition [6| . This condition relates the symmetry potential to the isovector nuclear 
density (the neutron-excess density) via the isovector (non-spin- flip) part of the p-h interaction (I26|) : 

v{r) = 2i^"(n"(r) - nP{r)). (45) 

Here, the neutron and proton densities are calculated with taking into account nucleon pairing: 



1 

2 . 
A 



'« = i;^L(2jA + l)«M(0, (46) 



1 r 

^A#Ao 



noddir) = -^ J2 (2JA + l)vlxl{r) + [(2jao - 1)0 + l\xlir) [>• (47) 



The mean Coulomb field Uc{t) is also calculated selconsistently within the Hartree approximation via the proton 
density P5|) . (|47p . Thus, the above-described mean field is a partially selfconsistent one. Nevertheless, the isospin 
symmetry of the model Hamiltonian is slightly violated due to approximate description of the nucleon pairing (the 
use of the "diagonal" approximation, truncation of the s-p basis, different values of the G° and G^ pairing strengths). 
These approximations can affect the quality of the calculated Fermi excitations, but they seem to be of rather minor 
importance for calculations of the GT excitations considered in this work. 

C. Choice of model parameters and the calculation scheme 

Recently, a partially selfconsistent mean field has been adopted to describe experimental single-particle spectra 
in the doubly-closed-shell nuclei ''^Ca, ^^^Sn, ^ospj-, gj ^g ^ starting point, the fully phenomenological mean field 
of Ref. [ij] was used. The deduced strength parameters C/q, Uso, F^ = Cf {C — 300 MeV ■ fm?) and also the 
diffuseness parameter a reveal a weak A-dependence (the value rQ—1.27 fm is taken the same for all nuclei [iSJ). This 
point allows us to use for the nuclei under consideration the properly interpolated values of the model parameters 
mentioned above (Table 1). Then, the partially selfconsistent mean field P^ - (|T7)) and the corresponding s-p level 
scheme are calculated to get a basis for the BCS solution of the nucleon pairing problem accordingly to Eqs. (I5^ - (IM1) . 
This problem is solved consistently with calculation of the symmetry potential and mean Coulomb field by making 
use of the neutron and proton densities (l46l) , (|47)) and a set of the pairing strength parameters Gj^ . The latter are 
then fitted to reproduce the experimental pairing energies Epair |14| : 



^-- ^' \£,{N + l) + \s,{N)-\f'- - ^ 



Spa^r{N) = ±{ -£b{N + 2) - -£f,{N + 1) + -£b{N) - -£b{N - 1) + -£b{N - 2) )■ (48) 



with "-I-" or "-" related to N(even) or N{odd), respectively. As £i, in this expression (the corresponding experimental 
values are taken from Ref. [lal), we use the properly modified ground-state energies of Eqs. ([BS]) . (pG ]) . For the 
subsystems having Nc ± 2 nucleons {Nc is a magic number) we use a simplified expression for the pairing energy: 

^biN, ±l) + ^£b{N, ± 2) - ^£b{N, ±3) + ^ 



£pa^r{Nc ± 2) = --£b{N, ± 1) + -£b{N, ± 2) - -£b{N, ±3) + -£b{N, ± 4). (49) 



The parameters G" „ = Cg^ „, fitted as described above, are listed in Table 1. To reduce an artificial violation of the 
isospin symmetry, caused by the use of only discrete (quasidiscrete) s-p states in description of the pairing problem, 
we use the same number of these states Nd+qd for the proton and neutron subsystems. The Nd+qd values for the 
nuclei under consideration are also given in Table 1. The energy and the radial wave function of a quasibound state A 
are found using the representation of the radial Green function, which is valid inside the nucleus at the energies close 
to the quasibound state energy: 

9{\) {r, r ,e) = -^-— , (50) 

e — ex + Tri x 
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where T\ is the quasibound-state escape width. 

After the consistent solution of the pairing problem, we turn to evaluation of the GT* strength functions and 
amplitudes Mq^ within the present pn-cQRPA approach. To calculate these quantities in accordance with Eqs. 
p2| . mi), ([15]) and Eqs. dfT]), (dH), dH]), respectively one has to choose the strengths F^ = Cg' and G^ (or 
gpp = 2G^/(G5^ + Gp)) of the spin-flip parts of the p-h and p-p interactions. Using in evaluation of the GT~ strength 
function the realistic value gpp = 1, we find the strength g' to reproduce in calculations the experimental GTR energy. 
The fitted strength values are listed in Table 2. These values are weakly dependent on the variation of gpp in a vicinity 
of the point gpp — 1. Having chosen g', the amplitudes Mq!^, their components Mqj. and M^rp are calculated as 
fimctions of gpp for the pairs of nuclei ™Ge-Se, ^™Mo-Ru, -^^^Cd-Sn and ^'^°Te-Xe. Since the ^^^Sn nucleus has the 
filled proton shell, in evaluation of the amplitude Mq^j, for the pair ^^^Cd-Sn we identify the vacuum |0') with that for 
|0). The strengths gpp ^ 1, which allow us to reproduce in calculations the experimental Mq^ values, are also listed 
in Table 2. These strengths are further used to evaluate the 2i//3/3-deca.y GT reduced strength function S'o'rp^Ld) = 

uj~^SQrp ' {(jj)/M'qj, (normalized to unity) and the corresponding reduced running sum r^j,(uj) = J" S'Qj,{uj)duj (goes 
to unity at large w). Along with the above-mentioned decomposition these quantities show how the amplitudes Af^^ 
are formed from the partial contribution of 1^ states in the intermediate nucleus. Finally, the fitted g' and gpp values 
are used to calculate the GT^ strength functions. To compare with the experimental GT~ strength distribution it is 
convenient to consider the reduced running sum r^ (w) — J S^ji (uj)duj / S^rjij^, where Sqj^j^ is the GT^ strength of 

the GTR. We show also the (non-reduced) running sum Wq^{uj) — J" S^r^{u})du} for the GT+ strength in the final 
nucleus. All the calculations are performed without taking into account the quenching of the total GT strength. 

IV. CALCULATION RESULTS 

Within the pn-cQRPA approach described in the preceding Sections we calculate the 2z//3/3-decay amplitudes Mq^ 
and the corresponding GT* strength functions for the pairs of nuclei ^^Ge-Se, ^™Mo-Ru, ^^^Cd-Sn and ^'^"Te-Xe. The 
results are presented in Figs. 1.-4. Each figure contains: a) the amplitude Mq^ and its decomposition in two terms in 
accordance with Eq. (PT|) calculated as functions of the strength gpp-, b) the reduced 2i'/3/3-deca,y GT strength function 
S'Qrp[ijj) calculated at the gpp value fitted to reproduce in calculations the experimental amplitude M^^ (Table 2.); 
c) the reduced 2://3/3-decay running sum rQ^(w) calculated at the fitted gpp value; d) the calculated GT^ strength 

function S'^j;(a;) for the initial nucleus, e) and f) the running sums rQrf and Rqx for the initial and final nuclei, 
respectively. The results are compared with available experimental data. 

The calculated results shown in Figs, l.a), b), c) - 4. a), b), c) demonstrate how the amplitudes M^^ are formed for 
the nuclei under consideration. The general features are the almost linear dependence of the component Mqj, on the 
p-p interaction strength gpp and vanishing of these amplitudes at gpp k, 1. Thus, we confirm the conclusion of Ref. [5| 
(obtained within a simpler model) that the component Mqj^ is only determined by the p-p residual interaction and 
goes to zero at gpp = 1 where the spin-isospin SU(4)-symmetry is restored in the p-p sector of this model Hamiltonian 
^. As compared with Rcf. 5], in our calculations the components Mqj. (and, therefore, the amplitudes M^rp) are 
somewhat smaller since the strength of the mean-field spin-orbit term used in our calculations is smaller. Again we 
confirm the conclusion of Ref. Q that at gpp ~ 1 the component Mqj, essentially depends on the spin-orbit term 
only, which is the main source of SU(4)-symmetry breaking in nuclei. Details of the formation of the amplitudes Af^^ 
are dependent on the shell-structure of nuclei. As follows from Figs. 2.b), 2.c) and 3.b), 3.c), population of the 1+ 
ground state of the intermediate nuclei ^°°Tc and ^^^In dominates the formation of the amplitudes A/q^ for the pairs 
-'^™Mo-Ru and ^^^Cd-Sn, respectively. Within the used calculation scheme the proper GT^ and GT+ transitions are 
mainly to the single-quasiparticle gg ,3 — >■ gy/a and g^,^ — )■ gg,^ transitions, respectively. For the pairs ^^Ge-Se and 

^■^"Te-Xe the amplitudes Mq^ are formed by population of many 1+ states of the intermediate nuclei ^^As and ^"^^I, 
respectively (Figs, l.b), l.c) and 4.b), 4.c)). It means, in particular, that the intermediate states having a relatively 
large excitation energy and, therefore, a small B(GT+) value (like the GTR) can, nevertheless, play essential role in 
the formation of the amplitude Mq^. For this reason, the experimental studies aimed to deduce the corresponding 
B(GT='=) values for the intermediate 1+ states and then to reconstruct the amplitudes M^.^ seem to be sufficient in 
the case of the approximate single-state dominance. The proper examples are the studies performed for 1"^ states in 
^™Tc [r?| and "''in [3], but it might not be fully correct in the case of the study of ^'^As [i|. 

Turning to calculations of the GT^^-* strength functions, we start from the description of the GTR in the intermediate 



^ In the QRPA calculations using a realistic G-matrix residual interaction il9il this point of crossing zero is shifted to gpp < 1 due to the 
effect of SU(4)-symmetry breaking by the tensor forces in the residual interaction 
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TABLE I: The phenomenological mean field parameters, singlet p-ii and p-p interaction strengths, and the number of the bound 



and quasibound s-p states used in calculations. 



Nuclei 


f/o, MeV 


Uso, MeV-fm^ 


a, fm 


/' 


90,n 


9o,p 


Nb+qb 


''''Ge-'^^Se 


51.3 


34.4 


0.60 


1.17 


0.41 


0.32 


16 


ioomo-IOORu 


51.5 


34.2 


0.61 


1.13 


0.45 


0.41 


16 


"«Cd-"«Sn 


51.6 


34.1 


0.62 


1.06 


0.39 


0.33 


22 


i30Te_i30xe 


51.7 


34.0 


0.63 


1.09 


0.36 


0.36 


22 



TABLE II: The phenomenological triplet p-h and p-p interaction strengths. The experimental GTR energies from Refs. [la|-[la] 
and the amplitudes Mq^ from Ref. 



^ are also 


given. 








Nuclei 


E^,GTR, MeV 


9 


Mgt, MeV"' 


9pp 


'"Ge-'^Se 
ioomo-IooRu 
"«Cd-"«Sn 
i30,pg_i30xe 


11.13 [16] 
13.3 [17] 
10.04 [18] 
13.59 Lie] 


0.81 
0.92 
0.77 
0.88 


0.14 
0.24 
0.13 
0.03 


0.62 
0.91 
1.07 
0.99 



nuclei ^^As, ^""Tc, ^^^Sb, ^^°I. The strength functions ^^^(a;) are shown in Figs. l.d)-4.d) (similarly to the strength 
functions shown in Figs. l.b)-4.b) we add a small imaginary part to the s-p potential while calculating the s-p Green's 
function in order to make the presentation more visible). Using the experimental data of Refs. |16l - [l8| on charge- 
exchange reactions, we fit the p-h strength parameter g' to reproduce in calculations the experimental value of the 
centroid of the GTR energy (Table 2). For A=116 systems we have performed the fitting of the GTR in ^^^Sb in view 
of the absence of appropriate experimental data for ^^^In. However, the value B(GT'^^'')=0.26±0.02 has been deduced 
in Ref. [2] from the ^^^Cdfp,n)-reaction with population of the l"*" ground state in ^^^In. Within the interval Ex < 3 
MeV the calculated GT^^'-strength distribution in ^-'^^In exhibits one l"*" state, corresponding to the back-spin-flip 
transition 15^/3 ^ I59/2 ^^^° ^^^ ^^ ground state of "^In, with the value B(GT(-))=1.05 (for "6Sn-j.ii6s^j ^j^^jg 
transition is Pauli blocked). Bearing in mind that the low-energy part of the GT^^' strength distributions cannot be 
described in details within any pn-QRPA-based approach (for instance, in view of existence of the quasiparticle-phonon 
coupling), we compare with experimental data the running sums rQj.{uj), Rqj,(^uj). The calculated and experimental 
running sums r^j! shown in Figs. l.d)-4.d), at least qualitatively, agree. With the exception for the running sums 

Rqx for the pairs ^^Se-^^As (Fig. l.f)) and ^^^Sn-^^^In (Fig. 3.f)), a reasonable description of the data is found 
for the B{GT~^) values corresponding to the ground-state to ground-state 1+ transitions ^™Tc-^™Ru (Fig. 2.f) and 
116gj^_116jj^ (Fig. 3.f). 

V. CONCLUSIONS 

We have given a detailed formulation of a version of the proton-neutron continuum-QRPA approach that makes 
use of: 

(i) realistic zero-range interactions in the particle-hole and particle-particle channels; 
(ii) a modern version of the phenomenological isospin-self-consistent nuclear mean field; 

(iii) the full basis of the particle-hole excitations, and a rather large basis of the particle-particle ones (with inclusion 
of several quasi-bound s-p states). Within the approach we have studied "anatomy" of the Gamow- Teller 2v/3/3-decay 
amplitude for a number of nuclei. The concept of the broken spin-isospin SU(4)-symmetry in nuclei has been made 
use of in the study of the ^pp-sensitivity of M^^. A reasonable description of the Gamow- Teller strength functions in 
/3^-channels for the nuclei in question has been obtained within the approach. 

Sensitivity of the calculated observables to realistic variations of the mean field parameters is going to be studied 
elsewhere. 
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FIG. 1: a) the 2///3/?-decay GT amplitude Mq^ (full line) and their components Mqj, (dotted-point line), A/gy (dotted line) 
calculated as a function of the strength Qpp for the pair ^^Ge-Se; b) and c) the 2i^/3/3-decay reduced GT strength function and 
running sum, respectively, calculated at the strength Qpp from Table 2 for the same pair of nuclei; d) and e) the GT~ strength 
function and reduced running sum, respectively, calculated at the strengths g' and Qpp from Table 2 for the pair ^^Ge-As (full 
line), the experimental data (dashed line) are taken from Ref. [la]; f) the GT"*" running sum calculated at the strengths g' and 
gpp from Table 2 for the pair ^''Se-As (full line), the experimental data (dashed line) are taken from Ref. ^. 
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FIG. 2: a)-c) the same as in Figs. l.a)-l.c) but for the pair ^"''Mo-Ru; d) and e) the same as in Figs, l.d) and l.e) but for 
the pair ^"^Mo-Tc, the experimental data are taken from Ref. [131 ; f) the same as in Fig. l.f) but for the pair ^™Ru-Tc, the 
experimental data (noted by the circle) are taken from Ref. \y\ . 
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FIG. 3: a)-c) the same as in Figs. l.a)-l.c) but for the pair ^^"Cd-Sn; d) and e) the same as in Figs, l.d) and l.e) but for the 
pair ^^®Sn-Sb, the experimental data (related to the differential (^He,t)-reaction cross section at 0°) are taken from Ref. [10]; 
f) the same as in Fig. l.f) but for the pair ^^^Sn-In, the experimental data are taken from Ref. Q]. 
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FIG. 4: a)-c) the same as in Figs. l.a)-l.c) but for the pair ^'^"Te-Xe; d) and e) the same as in Figs, l.d) and l.e) but for the 



pair Te-I, the experimental data are taken from Ref. [l6|; f) the same as in Fig. l.f) but for the pair Xe-I 



